Misoriented bilayer graphene in a magnetic field: Optical transitions at commensurate 

angles 
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Misoriented bilayer graphene with commensurate angles shows unique magneto-optical properties. 
The optical absorption spectra of such a system strongly depend on the angle of rotation. For 
a general commensurate twist angle the absorption spectra has a simple single-peak structure. 
However, our studies indicate that there are special angles at which the absorption spectra of 
the rotated bilayer exhibit well developed multi-peak structures. These angles correspond to even 
symmetry of the rotated graphene with respect to the sublattice exchange. Magneto-spectroscopy 
can therefore be a potentially useful scheme to determine the twist angles. 



Graphene, mechanically exfoliated from graphite 
display truly remarkable electronic properties Q, and 
holds an immense potential to become a key ingradient 
for a new generation of electronic devices. The dynamics 
of electrons in a single sheet of graphene, a hexagonal 
honeycombed lattice of carbon atoms is that of massless 
Dirac fermions with linear dispersion, chiral eigenstates, 
valley degeneracy, and unusual Landau levels in an ex- 
ternal magnetic field Bilayer graphene, on the other 
hand, show quadratic dispersion and the charge car- 
riers there are characterized as massive chiral fermions. 
Interestingly, epitaxial graphene which is thermally 
grown on the C face of the SiC substrate, as well as 
graphene grown by chemical vapor deposition (CVD) 6], 
are multilayer films and yet, quite surprisingly display be- 
havior similar to that of a single layer graphene • These 
systems are known to have a high degree of rotational 
misalignments j8J . Theoretical studies of turbostratic bi- 
layer graphene |9l-[l^ have indicated that in this case the 
interlayer coupling is suppressed and the systems can be 
roughly considered as two decoupled layers of graphene. 
At the same time due to the modulated nature [ll| of 
the interlayer transfer integral, these systems show quite 
rich low-energy physics, which strongly dep ends on the 
nature of commensurate stacking faults [i2[ . 

In a rotated bilayer graphene the rotational stacking 
fault is determined by an angle 9 of rotation of one layer 
relative to the other [see Fig. 1(a)]. Each layer consists 
of two sublattices, A and B, and is characterized by two 
primitive translational lattice vectors: a = a(l,0) and 
b = a(-l/2, V3/2), where a = 0.246 nm is the lattice 
constant. The commensurate rotation in a twisted bi- 
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FIG. 1: Misoriented graphene bilayer with angle of rotation 
9, shown schematically in real space (a) and reciprocal space 
(b). The black solid dots and red open dots are the (a) atomic 
positions in real space, and (b) K and K' points in the recip- 
rocal space in different layers. 



For the even commensurate stacking fault both A and B 
sublattice sites of the two layers are coincident at some 
point, while for the odd stacking fault only A sublat- 
tice sites of two layers are coincident. The AA-stacking 
and Bernal stacking correspond to even and odd stacking 
faults with angles 9 = and 60°, respectively. 

The reciprocal lattice of a graphene layer consist of K 

where 
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and K' sets of points: K + G^,., K' + G, 

k ~ raGi -f kG2, m and k are integers, G^ = 
27r/a(l, l/-\/3) and G2 = 27r/a(0, 2/-\/3) are primitive 
reciprocal lattice vectors, and K = 27r/a(l/3, l/VS), 
K' = 27r/a(2/3, 0). For a rotated layer the reciprocal 
lattice is rotated by an angle 9 around the origin [see 
Fig. 1(b)]. Then the even commensurate stacking fault 
corresponds to the twist angles at which the K points 
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of reciprocal lattice of two layers are coincident 
K + G^ j^ = K{9) + G^, (9) at some values of k, m, k' , 



layer is defined by the condition [Tl|, f„^j. = ma + kb = and m'. For oddjitacking fault the K and K' points are 



k'b', where a' and b' are given by the coincident, i.e., K + G 



Tm'k' = 

rotation of the primitive vectors a and b by 9. The an- 
gles corresponding to the commensurate stacking fault 
are determined from: cos9 = (Sq"^ —p^)/ {3q'^ +p'^), where 
q > p > are integers. There are two types of commensu- 
rate rotations that are distinguished by their symmetr" 
even or odd, with respect to the sublattice exchange 
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K'{9) + G^,A9) [1^ 
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In the case of the commensurate stacking fault the in- 
terlayer coupling becomes a periodically modulated func- 
tion of position, which results in the interlayer coupling 
determined only by the coincident points of the reciprocal 
lattice that are found from K + G^ j. = K{9) + G^, ,.,{9) 
or K + G„ = K'{9) + G.,„, ^,{9). Therefore the inter- 
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layer coupling is characterized by the Fourier transform 
of the interlayer potential, t(k), taken at points K+G^^^ j,: 

k ^ t{K + G,„ f.). As a result, the effective low en- 
ergy Hamiltonian of the twisted layer at commensurate 
condition takes the form [121 
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where tg = i,„ ^.e''', p„ ± = ± ip, is the electron 
momentum operator for layer a = 1,2, 10^ m/s. 

The phase angle (j> is determined from (p = 27r/3(m — fc), 
which can have possible values or ±2tt/3. For example, 
for 6 — 38.2° the phase angle is cj) — 27r/3. 

For the odd-twisted bilayer the Hamiltonian for all 
twist angles is identical to the Hamiltonian of a bilayer 
graphene with Bernal stacking, but with suppressed in- 
terlayer coupling. We found new and unique features for 
the even-twisted bilayer. While at = the Hamilto- 
nian of a twisted bilayer is similar to that of a bilayer 
with AA stacking, at cj) = ±27r/3 the Hamiltonian of a 
twisted graphene becomes unique. We show below that 
this Hamiltonian exhibits interesting features in positions 
of the Landau levels (LLs) and in magneto-spectroscopy 
of the rotated bilayer. 

Introducing a magnetic field in the Hamiltonians ([1])- 
(m by replacing the momentum operator p with p+eA/c, 
where A is the vector potential, we obtain the LL ener- 
gies of a twisted bilayer. For odd and even bilayers the 
interlayer tunneling introduces different coupling struc- 
tures. For the odd bilayer the LLs of individual layers 
belonging to different LL indices are coupled, while for 
the even bilayer the LLs of same indices are coupled. 

For the odd bilayer the LL energy spectrum is similar 
to that of a bilayer graphene with Bernal stacking 
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where n = 0, 1, 2, 



±1 corresponds to the con- 



duction and valence bands, respectively, and Si = ±1 
determines the splitting of levels due to interlayer cou- 
pling. Here Eq = Hvp/iQ and Iq — \J ch/eB. For each n 
there are two particle-like and two hole-like LLs. 

For the even-twisted bilayer the interlayer tunneling 
couples the LLs of the two layers with same LL indices, 
resulting in a splitting of the originally degenerate levels. 



FIG. 2: A few lowest LLs of three different types of misori- 
ented bilayer versus the magnetic field: (a) the odd-twisted 
bilayer; (b) the even-twisted bilayer with phase angle </!> = 0; 
(c) the even-twisted bilayer with phase angle (j> — 2tt/3. The 
interlayer coupling is set to 50 meV. Red and blue lines in (b) 
and (c) correspond to the LL index n = 1, 2 respectively. 



This splitting depends on (p. The energy spectrum of a 
twisted graphene in this case is 
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'2n+-^ + 2siV2n-^cos((/)/2), (3) 



where n — 0,1,2,..., Sq — ±1, Si — ±1. li cj) — 
then e — ±\/2neQ ± t„j j., which is a simple splitting of 
degenerate LLs of the two graphene layers. This energy 
spectrum is similar to that of A A stacked bilayer. If 

<j> = ±27r/3 then e = ±^2nel + 1^^^ ± v^e^t^ ,^. Here 

the energy splitting of degenerate LLs of two layers is less 
than the corresponding splitting for = 0. 

In Fig. 2, the magnetic field dependence of LLs for dif- 
ferent types of twisted bilayer are shown for interlayer 
coupling strength of 50 meV. For the odd graphene bi- 
layer [Fig. 2(a)] the energy spectra is similar to that of a 
bilayer with Bernal stacking. There is one level with ex- 
actly zero energy, which does not show any magnetic field 
dependence. For the even-twisted bilayer there are two 
levels without any magnetic field dispersion [Fig. 2(b)]. 
These levels correspond to the LL index n = where 
the energy of the LLs for all values of ^ is £ = ±t„j j,. 
For n > the LL spectrum depends on the value of (f). 
The LL spectrum of even-twisted graphene with (j) — 
[Fig. 2(b)] is similar to that of a bilayer with AA stacking. 
The LL spectrum of a graphene bilayer with = ±27r/3 
[Fig. 2(c)] shows an unique feature: there is a clear gap- 
like behavior in the spectrum with no LLs within the 
interval {^/3/2)t^^ > e > — (\/3/2)t^ ^ for all values of 
the magnetic field. 

The LL wavefunctions of the twisted bilayer also de- 
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FIG. 3: Intensity of optical transitions (absorption) between 
the valence (sq ; = —1) and conduction (sq j — 1) bands 
shown for even- twisted bilayer with phase angle (j) — 27r/3 
as a function of the LL index. The magnetic field is 0.5 Tesla 
for (a) and (b), and 2 Tesla for (c) and (d). The optical tran- 
sitions are between the initial states with Si^i = 1 [(a) and (c)] 
or i = — 1 [(b) and (d)] and the final states with f = 1 
(dots) or Si f = —1 (triangles). Inset: schematic illustration 
of optical transitions from the initial valence band states with 
LL index n and different values of Si to the final conduction 
band states with LL index n -I- 1 and different values of s^. 



the even-twisted bilayer is the same as in an isolated 
graphene layer [l^ . Therefore, as far as the FQHE is con- 
cerned, the even-twisted bilayer can be considered as two 
decoupled graphene layers for any twist angle • Since 
the odd-twisted bilayer is similar to the bilayer graphene 
with Bernal stacking, the FQHE in an odd bilayer is sim- 
ilar to a bilayer graphene with Bernal stacking [l5l |. 

The phase /3 in Eq. (|4]) depends on both the twist an- 
gle (through (/)) and the LL index n. Although this phase 
does not influence the interaction properties within a sin- 
gle LL, importantly, it can modify the optical transitions 
between different LLs. For light polarized along the x 
axis, the operator M of optical transitions is 
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where cr^ is the Pauli matrix. Then the selection rule for 
the optical transitions between the initial state i, and the 
final state f , is the same as for a single graphene layer, i.e., 
?if = rt; -I- ±1. We consider only the optical transitions 



to higher excited state, i.e., 



-I- 1, as those have 



higher frequencies and are perhaps easier to observe. The 
intensity of the corresponding optical transitions is 



pend on the type of commensurate stacking fault. Below 
we consider in detail only the case of even bilayer, which 
shows new and interesting features. The wavefunctions 
for the even-twisted bilayer corresponding to the LLs ([3]) 
are 
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where V'„ is the wavefunction of the conventional LL with 
index n, /„ = 1,C„ = 1/2 if n > and /„ = 0,C„ — 
1/^/2 if n = 0. The phase /3 is defined as 



/3 — arcsin 
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which is nonzero only for nonzero values of (p. This il- 
lustrates the sensitivity of the LL functions of a twisted 
bilayer to the twist angle. 

Armed with the wavefunctions ^ we are now ready to 
evaluate the strength of the electron-electron interaction 
within a single LL. This interaction is characterized by 
the Haldane pseudopotentials 16] and is responsible for 
formation of the fraction quantum Hall effect (FQHE) 
states [l3|. Since the interlayer coupling in Eq. ^ af- 
fects only the phases of the wavefunction components, 
the pseudopotentials for the wavefunctions ^ are de- 
termined only by the LL index n and does not depend 
on the interlayer coupling. The pseudopotentials for the 
even-twisted bilayer are identical to those of individual 
graphene layers. Hence the strength of the FQHE in 
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For the even- twisted bilayer with = the phases /3; and 
/3f are zero. Then the intraband optical transitions (the 
same sign of Sq) are allowed only between the states with 
the same sign of Si, while the interband optical transi- 
tions are allowed between the states with opposite signs 
(4) of s^. All transitions have the same intensity. Therefore 
for a given LL index ri; the optical absorption consists of 
a single line. 

A different situation occurs for the even-twisted bilayer 
with nonzero values of (j). In this case both /3j and /3f are 
nonzero and depend on the LL index n. We consider be- 



low only the interband optical transitions, i.e., Sq i ~ "-'^ 
in the initial state (valence band) and Sg j. = 1 in the final 
state (conduction band). For each initial state there are 
two transitions to the final states with si^f = 1 and — 1. 
In Fig. 3 the relative intensity of these transitions are 
shown for different initial LLs and for different magnetic 
fields. For each LL index the initial and final states are 
characterized by the parameter s^. The four possible op- 
tical transitions at a given LL index are shown schemat- 
ically as inset in Fig. 3(c). In Fig. 3(a,c), the transi- 
tions from the states with j = 1 are shown, while in 
Fig. 3(b,d) the initial state has s-^ ; = —1. The triangles 
and dots correspond to the final states with s-^ ^ = — 1 and 
1, respectively. The general tendency illustrated in Fig. 3 
is the existence of strong transitions to both final states. 
At a small magnetic field and small LL index, these tran- 
sitions have comparable intensities [see Fig. 3(a,b) where 
the magnetic field is 0.5 Tesla]. With increasing LL index 
the transition to one of the states is suppressed and the 
system becomes similar to the case oi cj) = 0. The regions 
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FIG. 4: Intensity of optical transitions (absorption) be- 
tween the valence (so.i = ~1) and the conduction (sq f = 1) 
bands is shown for the even-twisted bilayer with a phase angle 
= 27r/3 as a function of the transition energy for different 
LL index rii of the initial state. The magnetic field is 0.5 Tesla 
(left panel) and 2 Tesla (right panel). The black lines corre- 
spond to transitions from the initial state with = 1 only, 
while the red lines correspond to transitions from all initial 
states (si = 1 and —1) with a given LL index n. For compar- 
ison, the optical absorption spectra of odd bilayer (blue line) 
and even bilayer with (j) = (green line) are shown in panel 
(b) for Hi — 1. 

of magnetic fields and LL indices for which the optical 
transitions have comparable intensities are determined 
by the strength of inter-layer transfer integral, tm,k- The 
strongest coupHng should be expected for small values of 
m and k, e.g., for the twist angle 38.2°. 

The unique features of optical transitions in the even- 
twisted bilayer with are clearly visible in the optical 
absorption spectra. In Fig. 4, we show the absorption 
spectra from a given LL of the valence band of a twisted 
bilayer with interlayer coupling of 50 meV. At a given 



LL index the absorption spectra shown by black solid 
lines correspond to transition from the ground state with 
s-^ ; = 1 [see the inset in Fig. 3(c)], i.e., only the state with 
; = 1 is occupied. The absorption spectra from all the 
initial states (s^ j = 1 and —1) with a given LL index are 
shown by red lines. For the LL index n = there is only 
one initial state in the valence band. A 2 meV broadening 
of the optical lines is introduced in Fig. 4. The absorption 
spectra clearly show a multi-peak structure which is more 
pronounced at a small magnetic field [0.5 T in Fig. 4(a-c)] 
and at small LL index. With increasing magnetic field 
[see Fig. 4(d-f)] or increasing LL index only one peak 
in the absorption spectra survives, which is consistent 
with the results shown in Fig. 3. Such a behavior of the 
absorption spectra of even-twisted bilayer with 7^ is 
totally different from the odd-twisted graphene or even- 
twisted with = 0. It is easy to show that for the odd- 
twisted graphene, which is similar to bilayer graphene 
with Bernal stacking, only one strong optical transition 
exists for each LL. Therefore for both odd bilayer and 
even bilayer with (f) — the optical spectra for each LL 
consist of a single line [Fig. 4(b), ni = 1]. Figure 4 clearly 
shows the fingerprints of phase angle in the magneto- 
optics of commensurate twisted graphene. 

In conclusion, magneto-optical properties of twisted 
graphene bilayer show strong dependence on the twist 
angle. At twist angles corresponding to odd bilayer and 
even bilayer with zero phase angle, = 0, the absorption 
spectra from a given LL consist of a single line, while the 
optical spectrum of the even bilayer with nonzero angle (j) 
has well-developed multi-peak structure which can per- 
haps be observed experimentally. Such dependence of 
the optical spectra on the twist angle is visible at low 
LL index, n < 10, and weak magnetic field, B < 5 Tesla. 
The strongest multi-peak structure should be observed at 
large inter-layer coupling, e.g. for the twist angle 38.2". 
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